• Review • 



SCIENCE CHINA 

Physics, Mechanics & Astronomy 



December 2012 Vol. 55 No. 12: 1-14 
doi: 10. 1007/sl 1433-012-4929-9 



Three-dimensional topological insulators: A review on host materials 

FENG Wanxiang and YAO Yugui* 

School of Physics, Beijing Institute of Technology, Beijing 100081, China 
Received August 30, 2012; accepted October 22, 2012 



In recent years, three-dimensional topological insulators (3DTI) as a novel state of quantum matter have become a hot topic in the 
fields of condensed matter physics and materials sciences. To fulfill many spectacularly novel quantum phenomena predicted in 
3DTI, real host materials are of crucial importance. In this review paper, we first introduce general methods of searching for new 
3DTI based on the density-functional theory. Then, we review the recent progress on materials realization of 3DTI including simple 
elements, binary compounds, ternary compounds, and quaternary compounds. In these potential host materials, some of them have 
already been confirmed by experiments while the others are not yet. The 3DTI discussed here does not contain the materials with 
strong electron-electron correlation. Lastly, we give a brief summary and some outlooks in further studies. 
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1 Introduction 

In the past few decades, topological quantum states charac- 
terized by specially topological orders [1,2] have become one 
of the most critical physical phenomena in condensed mat- 
ter physics. In general concept, topological insulator con- 
tains two classes of topological quantum states depending on 
whether the system has time-reversal symmetry (TRS) or not. 
In the case of two-dimension, the first is the quantum Hall 
(QH) state [3], in which the TRS is explicit broken due to 
extemal magnetic field or intemal magnetization. This topo- 
logically nontrivial state can be classified by the first Chern 
number proposed by Thouless et al. [2], directly connected to 
the quantized Hall conductivity. The second is the quantum 
spin Hall (QSH) state [4-6 J, which is protected by TRS and 
essentially resulted from spin-orbit coupling (SOC). Similar 
to the role of the Chern number, Z2 topological invariant in 
TRS protected systems was introduced by Kane and Mele [7], 
that is, Z2=l (0) represents topologically nontrivial (trivial) 
state. The topological insulator with the common feature, 
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having gaped insulating states in the bulk but gapless surface 
states on the edge, fundamentally distinguishes itself from the 
ordinary insulating state originally defined by Kohn [8]. 

The QSH state, i.e., two-dimensional topological insulator, 
has attracted research focus because of its promising applica- 
tions in spintronic devices. The existence of QSH state in 
graphene was first proposed by Kane and Mele [5,6], how- 
ever, subsequent works [9-11] showed that the band-gap in- 
duced by the SOC is too small (~ 10"'' eV) and QSH state in 
graphene can not be observed in current experimental condi- 
tions. A more realistic material is the HgTe/CdTe quantum 
well, theoretically predicted by Bernevig et al. [12] and then 
experimentally confirmed by Konig et al. [13]. Soon after the 
QSH state was discovered, topological quantum states under 
TRS were also extended in three -dimension [14-16], i.e., the 
three-dimensional topological insulators (3DTI). In 3D, there 
are four Z2 topological invariants, vq; (V1V2V3), used to clas- 
sify the topology of an insulator [17]. A nonzero vq indicates 
that the system is a strong topological insulator (STI). When 
vo=0, the systems are further classified according to vi, V2, 
and V3. The systems with yi,2,or3=0 (not all of them) are 
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called weak topological insulators (WTI), while 0;(000) is a 
normal insulator (NI). A STI can not be directly connected to 
a WTI or a NI by any adiabatically continuous transforma- 
tion of band structure. In the WTI or NI, the surface states 
have an even number of Dirac points and can be violated by 
structural disorders or ferromagnetic impurities. Conversely, 
in the STI, the surface states have an odd number of Dirac 
points and are topologically protected. 

The 3DTI is fascinating mainly because of its exotic sur- 
face states [18-21], which are very robust and against from 
the extrinsic perturbations, such as structural disorders and 
nonmagnetic impurities. Combining the topological surface 
state with superconducting state or ferromagnetic state, the 
3DTI has shown great application potentials in quantum com- 
puting and spintronics fields [22-23]. So far many 3DTI 
host materials have been theoretically predicted and/or ex- 
perimentally confirmed [see Table [T]. However, most of cur- 
rent 3DTI host materials have various disadvantages as ob- 
served in experimental observations, for example the large 
concentration of bulk carrier, which hinders further investi- 
gation along the direction of practical applications. There- 
fore, searching for new better host materials with simpler 
electronic structures and simpler synthetic conditions is a key 
step in this field. 

In this review paper, we focus on the 3DTI host mate- 
rials. We introduce general methods of searching for new 
3DTI based on the density-functional theory. We also discuss 
many theoretically proposed 3DTI materials including sim- 
ple elements, binary compounds, pseudo-binary compounds, 
ternary compounds, and quaternary compounds. In these po- 
tential host materials, some of them have already been con- 
firmed by experiment while the others are not yet. The 3DTI 
discussed here does not contain the materials with strong 
electron-electron correlation. 

2 How to search for 3D-topological insulators 

As a first step towards various appUcations of 3DTI, mate- 
rial realization is of crucial importance. During the search 
of 3DTI in real materials, first-principles calculations guided 
by the topological band theory [18,19] have played a pivotal 
role. We briefly introduce four widely used methods to search 
3DTI. In practice, two or more among these methods are usu- 
ally combined together to determine the band topology of a 
crystal. 

2.1 Surface state electronic structure 

As a fundamental feature in 3DTI, there is gapless surface 
state inside the bulk band-gap [18-21]. The gapless surface 
state has the form of linear Dirac cone, which holds massless 
relativistic particles described by the Dirac equation [24,25]. 
Furthermore, the spin and momentum of surface Dirac elec- 
trons are locked together in the sense that the electrons with 
opposite spin must propagate in opposite direction along the 
boundary due to TRS. Hence, one can determine whether an 



insulator is a 3DTI or not by calculating the surface state elec- 
tronic structure. An odd number of surface states crossing the 
Fermi level indicates a 3DTI, otherwise it is a normal insu- 
lator At the experimental aspect, the surface state electronic 
structure can be directly measured by the angle-resolved pho- 
toemission spectroscopy (ARPES). The comparison between 
theoretical calculations and experimental ARPES measure- 
ments becomes possible. 

Although the bulk-boundary correspondence as a basic 
routine can be used to find the topologically nontrivial mate- 
rials, it is not eff'ective in first-principles calculations. Physi- 
cally speaking, the surface band dispersion is sensitive to the 
details of the surface, such as the terminations and orienta- 
tions [26-28]. Furthermore, some topologically trivial states, 
for example the dangling bonds on the surface of group-IV 
or III-V semiconductors, may manifest themselves in the sur- 
face band structure. The surface reconstruction or hydro- 
gen passivation can partly destroy these dangling bonds, but 
the coexistence of topologically trivial and nontrivial surface 
states certainly complicates the identification of the topo- 
logical order From the computational point of view, first- 
principles calculation of surface state electronic structure is 
an inefficient method because it generally requires a large 
amount of computational resources. 



2.2 Adiabatic continuity band transformation 

The argument of adiabatic continuity was once used to iden- 
tify the quantum spin Hall phase in HgTe/CdTe quantum well 
[12] and silicene [29], but can also be introduced into the 
3DTI. The basic idea is that the Hamiltonian of a system can 
be adiabatically connected to that of another system. For the 
former and latter systems, the topological invariants can only 
change when the bulk band-gap closed. Therefore, if the band 
structure of a material can be smoothly transformed to that 
of another one without closing the band-gap, then they must 
share the same topological classification. In contrast, if the 
band-gap closes during the band transformation, the topolog- 
ical phase transition may occur. 

In practice, the adiabatic continuity transformation of band 
structure can be artificially controlled, for example, by al- 
tering the composition x in solid alloy Bii_vSb^ [17,30,31] 
or Pbi_^SnfTe [17]. By starting x from zero, these systems 
undergo the topological phase transition from trivial state to 
nontrivial state. If we continue increasing x for Pbi-j^Sn^Te, 
a further transition from nontrivial state to trivial state will 
happen. Tuning other parameters, such as atomic number 
and strength of SOC, are also widely used in literatures. 
The obvious advantage of this method is its apparent ease. 
One can predict new 3DTI candidates based on some well- 
known topologically nontrivial or trivial materials. While 
many intermediate phases are required along the transforma- 
tional path connecting the initial state and the final one, which 
makes this method computationaUy demanding. 
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2.3 Band inversion picture 

Sometimes, topologically nontrivial materials can be well 
recognized at the first glance by the so-called band inver- 
sion picture [12]. In most of the common semiconductors, 
the valence band-edges are formed by the p-orbits of elec- 
trons, whereas the conduction band-edges are formed by the 
s-orbits of electrons. This situation belongs to normal band 
order, i.e., topologically trivial phase. In other cases, the rel- 
ativistic effect from heavy elements can be so large that the s- 
orbits is pushed below the p-orbits, that is, the inverted band 
order appears. A typical example is HgTe in which the effec- 
tive positive charge of Hg core is increased due to the partly 
delocalization of its d-orbits. This causes the i-orbits to be 
more attracted by the Hg core, and consequently, the i-orbits 
are pushed below the p-orbits. i.e., the inverted band struc- 
ture forms [32]. 

The band inversion is a strong indication that a material 
falls into topologically nontrivial phase. The first application 
is the prediction of quantum spin Hall effect in HgTe/CdTe 
quantum well proposed by Berne vig et al. [12]. The quan- 
tum well structure is that the HgTe with inverted band order 
is sandwiched by the CdTe with normal band order. The topo- 
logical property of entire quantum well is determined by the 
thickness d of HgTe layer. The critical thickness dc is pre- 
dicted to be about 6.5 nm [12]. For a thin well when d < dc, 
the CdTe has the dominant effect and the entire quantum well 
is topologically trivial with normal band order, while for a 
thick well when d > dc, the HgTe has a critical role role 
and the entire quantum well is topologically nontrivial with 
inverted band order. 

Similar to HgTe and CdTe, the band topology of other cu- 
bic semiconductors with zinc-blende-like structure can also 
be determined by the band inversion picture. To do this, one 
need to define a useful physical quantity, namely, the band 
inversion strength (BIS). Specially at F point, the BIS is de- 
fined as the energy differences between Fg state (formed by 
i-orbits) and Fg state (formed by p-orbits) [33,34] i.e., 

AE = Er, - Er, . (1) 

In general, negative AE typically indicates that the materials 
are in topologically nontrivial phase, while those with posi- 
tive AE are in topologically trivial phase. 

Moreover, the conventional jp-orbits semiconductors, the 
band inversion has also been found in pp-orbits and <i/-orbits 
semiconductors. In Bi2Se3 family, both the valence and con- 
duction band-edges are formed by p-orbits but with oppo- 
site parities. It was predicted that Bi2Se3 family is 3DTI due 
to the inverted band order between \Pit^ and \P2^^ orbits 
[35]. Here, the + (-) means the parity of the corresponding 
Bloch wavefuntions. The band inversion can be interpreted 
that at F point these two p-orbits with opposite parities are 
exchanged when SOC turns on. Other examples are actinide 
compounds AmX (X=N, P, As, Sb, and Bi) and PuF (y=Se 
and Te), which were recently discovered by Zhang et al. [36] 



as a new class of 3DTI driven by strong electron-electron in- 
teractions. In these compounds, the Fg state of 6(i-orbit lo- 
cates below the F^ state of 5 /-orbit at X point, and as a result, 
the inverted band order occurs. 

Although the band inversion at some high-symmetry 
points is a convenient way to find topologically nontrivial ma- 
terials, for example HgTe/CdTe quantum well, BiiSea family, 
and actinide compounds, it should be used with more care be- 
cause the band topology is a global property within the entire 
BriUouin zone and is not only Umited to some high-symmetry 
points. 

2.4 Z2 topological invariants 

The most general and direct method of searching for 3DTI is 
to calculate the Z2 topological invariants from the bulk band 
structure [14-16]. There exist two situations depending on 
whether the system has spatial inversion symmetry or not. If a 
system has inversion symmetry, the calculation of Z2 topolog- 
ical invariants can be well simplified based on the parity cri- 
terion developed by Fu and Kane [17]. Conversely, if a sys- 
tem does not have inversion synmietry, Fukui and Hatsugai 
[37] have developed an effective algorithm, which requires 
Bloch functions (BFs) on a dense two-dimensional grid to 
compute the Z2 topological invariants. The implementation 
of these two methods within full-potential linearized aug- 
mented plane-wave (FP-LAPW) formalism has been demon- 
strated [38]. In the following, we briefly introduce these two 
methods. 

In the systems with inversion symmetry, the parity anal- 
ysis only requires the knowledge of BFs on the time- 
reversal invariant momenta (TRIM) in the Brillouin zone 
(BZ) [17]. In 3D system, there are eight TRIMs, r,=(„j„2„3) = 
2 (n\Gi + ^262 -I- n^Gi), where Gj are primitive reciprocal- 
lattice vectors with rij = or 1. The Z2 invariants are deter- 
mined by the quantities 

<J; = n^2m(r,). (2) 

m=l 

Here, ^2m (Xi) = <*l'2m,r, 1^1 ^imXi) is the eigenvalue of parity 
operator P at the 2m-th occupied band and TRIMs F,. The 
^2m(r,) is equal to 1 (-1), corresponding to even (odd) par- 
ity of the BFs. The sum is over all of the occupied bands 
with only even band index due to the Kramers degeneracy at 
TRIMs. In 3D system, there are four independent invariants 
vo; (V1V2V3), given by [17] 

8 

!=1 

(-1)"'= n ^M«,«2«3). (4) 

nt=l,ny*t=0,l 

where vq is independent of the choice of primitive reciprocal- 
lattice vectors Gj while vi, V2, and V3 are not. The combina- 
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tion of these four independent invariants vo; (vi V2V3) clearly 
distinguish three classes of states: STI, WTI, and NI [17]. 
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Figure 1 Schematic drawing of lattice mesh in a 2D-torus. Under the time- 
reversal constraint, only half of 2D-torus is needed, which is denoted by 
shaded region. The thick lines indicate the boundary of S* , i.e., dS^ , and 
the open arrows denote their directions. All fc-points are divided into three 
classes: B*, S;, and S°, which are represented by small solid, small open 
and large shaded circles, respectively. Adapted from Ref. [38]. 

In the systems without inversion symmetry, an effective al- 
gorithm on a dense two-dimensional grid to compute the Z2 
invariants has been proposed by Fukui and Hatsugai within 
a tight-binding framework [37]. It was shown that under the 
time-reversal constraint the Z2 invariants can be written in 
terms of the Berry gauge potential and Berry curvature asso- 
ciated with the BFs [39], 



Z2 = 



27T 



(D dk-Jlik)- f cfikTik) 

JdS* J'B* 



mod 2, (5) 



where J?l (k) and T (k) are the Berry connection and Berry 
curvature, respectively. 



and 



Jl(k) = /^<M„ (k) I ViM„ (k)) 



r{k)^VtxM(k) I, 



(6) 



(7) 



The and dS^ indicate half of 2D-torus and its boundary, 
respectively (see Fig. [TJ. By using the periodic gauge [40,41] 



\u„{k + Gj}) = e-'''r'\u„(k)}. 
and two time-reversal constraints [37,39] 



\u„(-k)) = &\u„(k)), ke 



(8) 



(9) 



|m2« (-^)> = \u2n-i (k)} , -A: and A: e . 



(10) 



one can obtain the periodic part of BFs \u„ {k)) at every k- 
point of the 2D-torus. Then, the Berry connection [Eq. ([6])] 
and Berry curvature [Eq. (jT])] can be calculated by the finite 
element expressions [37]. After that, the Z2 invariants can 
be obtained by inserting Eqs. (|6]l and (|7]i into Eq. Q. In 
3D system, there are six 2D-tori T(Zo), T{Zi), T(Xo), T(Xi), 
T(Yo), and T{Yi), supporting six Z2 invariants zo, Z\, xo, xi, 
yo, and yi. However, out of these six possible Z2 invariant 
only four of them are independent due to the constraint xq + 
x\ - yQ+y\ = Zo + -Zi (mod 2). The Z2 invariants are denoted 
in another way with vq = (zo + Zi) mod 2, v\ - x\, - y\ 
and V3 - zi, i.e., vo;(viV2V3) [14-16]. By using topological 
invariants vo; (V1V2V3), one can determine the band topology 
for a given material. 

Comparing with the surface state electronic struc- 
ture (Sec. |2.1| l, adiabatic continuity band transformation 



(Sec. 2.2 1, and band inversion picture (Sec. 2.3 1, the calcu- 



and 



lation of the Z2 topological invariants is a direct evidence for 
determining topologically nontrivial state. Furthermore, from 
the viewpoint of first-principles calculation, it is the most 
computationally efficient. Our implementation of the calcu- 
lation of Z2 topological invariants within FP-LAPW formal- 
ism has been recently applied to predict ternary half-Heusler 
[33,34] and chalcopyrite [42] 3DTI and quantum spin Hall ef- 
fect in silicene thinfilm [29]. Since the critical factor is to cal- 
culate the eigenvalues of the parity operator and time-reversal 
operator in the systems with and without inversion symme- 
try respectively, it can be implemented in any other first- 
principles methods, such as the pseudopotential planewave 
and the linear muffin-tin orbital. There appears another ap- 
proach to compute the Z2 invariants by employing the charge 
center of Wannier functions [43,44]. 

3 3DTI host materials 

We review the 3DTI host materials in the sequence of sim- 
ple elements, binary compounds, pseudo-binary compounds, 
ternary compounds, and quaternary compounds. In each part, 
only one or several typical materials are discussed. Table [T] 
includes a list of all 3DTI host materials. In these poten- 
tial host materials, some of them have been confirmed by ex- 
periments while the others are not yet. The readers can find 
the interesting materials in Table[T]and refer to the literatures 
therein. It also should be noted that the 3DTI discussed here 
does not contain the materials with strong electron-electron 
correlation. 

3.1 3DTI in simple elements 

3.1.1 a-Sn 

The group-IV element a-Sn (grey tin) crystallizes in the dia- 
mond structure with space group Fd3m (No. 227), as shown 
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in Fig. |2|a). The full-relativistic band structure in Fig. |2|c) 
clearly shows that at the static lattice constant or-Sn is a 
zero band-gap semiconductor with the feature of band inver- 
sion. The valence band maximum (VBM) and the conduction 
band minimum (CBM) are degenerated at F point and the 
Fermi level crosses the fourfold-degenerated F^ state (from 
/7-orbits). The i-orbits-like F^ state is situated below the Fg 
state, forming the inverted band structure [48,49]. Hence, the 
band structure of a-Sn is qualitatively different from those of 
the other group IV elements C, Si, and Ge [49]. 

Fu et al. [17] proposed that a uniaxial strain can tune a-Sn 
into a 3DTI. The basic idea is that after applying a uniaxial 
strain the fourfold-degenerated Fg state is lifted and a band- 
gap opens around F point. Consequently, the topologically 
nontrivial state forms because the band inversion does not 
change. Fig. |2|d) shows the band structure of a-Sn under 
a uniaxial strain c/a=l-i-3% but with constant volume. One 
can see that the band-gap opens while the band inversion pre- 
serves. To further confirm the argument of band inversion, 
the Z2 topological invariant is calculated. The result 1;(000) 
indicates that the strained a-Sn is indeed a STI. Although the 
theoretical picture is very clear, there is not any experimental 
work on a-Sn currently. 



(c) Static (d) strained 




strain c/a=1+3% 



WLrx PNrx 



Figure 2 Crystal structure of a-Sn at its static lattice constant (a) and 
under a uniaxial strain c/a=l+3% (b), respectively, (c) and (d) are the corre- 
sponding band structures with the i-orbits projections, i.e., the size of solid 
circles. The calculated Z2 topological invariant 1;(000) indicates a strong 
3DTI. Here, the full-relativistic band structures were performed using FP- 
LAPW method [45] with the MBJLDA exchange-correlation potential [46], 
implemented in wien2k package [47] . 

3.1.2 Sb 

The Sb (antimony) crystallizes in rhombohedral structure 
with space group R3m (No. 166), as shown in Fig. [3|a). 
The full-relativistic band structure in Fig. [3|c) shows Sb is a 
semimetals with small pockets of electron in the vicinity of 
the L point. There are also small pockets of hole around the 
low-symmetry H point (not shown). Although it is not an in- 
sulator, Z2 topological invariant can still be defined because 
the local band-gaps (negative indirect band-gap) exist at ev- 
ery ^-point throughout the entire BZ [17]. The result 1;(111) 



indicates that the Sb is a STI. As another group-V element, 
Bi (bismuth) has similar crystal and band structure compar- 
ing to Sb, as shown in Fig. |3ja) and[3]^d) respectively, but Bi 
is a NI with Z2 topological invariant 0;(000). The different 
topological classes between Sb and Bi were well interpreted 
by their different parities at L points [17]. 

Hsieh et al. [51,52] have measured the bulk and surface 
electronic structure of Sb. Because bulk Sb is a semimetal, it 
is difficult to separate its surface state from the projection of 
bulk state around the Fermi level. While by using advanced 
spin-resolved ARPES technique, Hsieh et al. [51] success- 
fully isolated the surface states of Sb(l 1 1) from its bulk states 
over the entire BZ. They also directly found the gapless and 
spin-splitting surface state, which characterizes the topologi- 
cally nontrivial feature in Sb [51]. 



(c) Sb 



(d) Bi 




Figure 3 Crystal structure (a) and the Brillouin zone (b) of Sb (Bi). (c) 
and (d) are the band structures of Sb and Bi, respectively. The Sb is a STI 
with Z2 topological invariant 1;(111), while Bi is a NI with 0;(000). Here, 
the full-relativistic band structures were performed using FP-LAPW method 
[45] with the PBE-GGA exchange-correlation potential [50], implemented 
in wien2k package [47]. 



3.2 3DTI in binary compounds 

3.2.1 Biy-xSbx alloy 

Both Sb and Bi are group-V semimetals with the rhombohe- 
dral crystal structures and similar lattice constants, therefore, 
they can easily form solid alloy Bii_iSbf [53]. By altering 
the concentration x of Sb substitution, Bii_jSbv undergoes a 
phase transition between semimetal and semiconductor In 
the range of 0.07 < x < 0.22, Bii_fSbv becomes semicon- 
ductor with the largest global band-gap of 0.03 eV at x=0.18. 
Based on a tight-binding model [54], Fu et al. [17,30] pre- 
dicted that Bii^jSbj; is the first realistic 3DTI material. Sub- 
sequently, this prediction was confirmed by first-principles 
calculation [31] and experimental observations [55-57]. 

Although both theories and experiments got a consistent 
conclusion that the topologically nontrivial phase exists in 
Bii_vSb^ alloy, the discrepancies about surface electronic 
structure still remain. The main concern is the surface band 
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configuration, including the numbers of surface bands and 
the crossing times between the surface bands and the Fermi 
level along the T-M line in surface BZ. All the experimen- 
tal works [55-57] reported three surface bands S2, and 
Ej) lying inside band-gap of bulk projection, but with differ- 
ent crossing times with the Fermi level, five times in Refs. 
[55,56] and three times in Ref. [57]. Conversely, both tight- 
binding [30] and first-principles [31] calculations can not re- 
produce the presence of the third surface band S'j. Teo et 
al. [30] and Zhang et al. [31] give three and five times of 
band crossing, respectively. Therefore, more experimental 
measurements and theoretical calculations are needed for re- 
vealing the microscopic physics in Bii_;cSb;c. 



(a) I \ (b) / \ (c) 



\ / \ MiT^\ / \ / 

\ J! .. ; ■-•—!! .\i \ / 



Figure 4 The band order of cubic diamond/zinc-blende semiconductors at 
the r point, (a) The normal band order, (b)-(c) The inverted band order can 
be obtained from the normal band order by either increasing the strength of 
SOC (b) or by decreasing the coupling potentials of .vi-orbits or pp-orbits 
(c). The band inversion strength is defined as AE = Ey^ - E^g, which has 
positive value in normal band order and negative value in inverted band or- 
der. Adapted from Ref [75]. 

3.2.2 Bi2Se^ family 

Since the surface electronic structure of Bii_fSbj is rather 
complicated and its bulk band-gap is too small, searching for 
new 3DTI host materials with simple surface electronic struc- 
ture and large band-gap becomes extremely important. For- 
tunately, BiaSes family compounds with larger band-gap and 
simper surface spectrum have been found as a second gen- 
eration of 3DTI materials [35,58-60], which support further 
study on various topologically protected phenomena at room 
temperature. 

Tetradymite semiconductors BiaSea, Bi2Te3, and Sb2Te3, 
have a rhombohedral crystal structure with space group R3m 
(No. 166). Although Sb2Se3 essentially has an orthorhom- 
bic crystal structure, a virtual rhombohedral crystal structure 
is used at here for comparing its topological property with 
tetradymite semiconductors. Zhang et al. [35] theoretically 
proposed that Bi2Se3, Bi2Te3, and Sb2Te3 are 3DTI while 
Sb2Se3 is not. More importantly, the band-gap of Bi2Se3 is 
as large as 0.3 eV, which is larger than the energy scale of 
room temperature. Almost at the same time, these family 
compounds were experimentally confirmed with the obser- 
vation of single Dirac cone on the surface for Bi2Se3 [58], 
Bi2Te3 [59,60], and Sb2Te3 [60]. Because of the large band- 
gap and simple surface electronic structure, Bi2Se3 family 



compounds have been extensively studied ranging from bulk 
tothinfilm [61-71]. 

3. 2. 3 HgTe and strained InSb 

HgTe is a zinc-blende compound with space group F43m 
(No. 216). It is a zero band-gap semiconductor due to large 
relativistic effect. The twofold-degenerated Fe state is located 
below the fourfold-degenerated Fg state, forming the inverted 
band order [see Fig. ^b)]. Fu et al. [17] theoretically pre- 
dicted that HgTe is a 3DTI under a uniaxial strain. The basic 
routine is similar to that of a-Sn, that is, the uniaxial strain is 
used to break the cubic symmetry and open a band-gap at F 
point without changing the band inversion. Recently, Briine 
et al. [72] successfully observed the Dirac-like topological 
surface state on strained HgTe by ARPES measurement. 

InSb is also a zinc -blende semiconductor, but in contrast 
to HgTe, it has a small band-gap of 0.235 eV [74]. The Fg 
state is located above the Fg state, forming the normal band 
order [see Fig. |4fa)]. Feng et al. [75] proposed that InSb can 
be turned into a 3DTI by a 2% ~ 3% biaxial lattice expan- 
sion. The generic guiding principle is that lattice expansion 
decreases the coupling potentials of ii-orbits or /?/?-orbits, 
which leads to the band inversion [see Fig. |4|c)]. Here, we 
use different modes of strain in HgTe and InSb. In the for- 
mer one, the uniaxial strain is just used to break the cubic 
symmetry. While in the latter one, the nonhydrostatic (2% ~ 
3% biaxial) lattice expansion is twofold because it not only 
changes normal band order to inverted band order but also 
breaks the cubic symmetry. 

3.3 3DTI in pseudo-binary compounds 

Although Bi2Se3 family as the second generation of 3DTI 
has attracted research focus because of its simple surface 
electronic structure, it hard to be applied in various real de- 
vices. The prominent difficulty is high bulk carrier concen- 
tration, which remarkably masks the contribution of surface 
carriers in the surface transport measurement [55,62,80-82]. 
The high bulk carriers concentration mainly originates from 
the native crystal defects [58-60]. Some experimental tech- 
niques including compensate doping in bulk [59,63] and tun- 
ing Fermi level by gate voltage [83-85] have been tested, but 
this issue is not perfectly solved. 

Recently, a promising platform for reducing the bulk car- 
rier concentration, Bi2Te2Se, has been proposed by Ren et 
al. [86] and Xiong et al. [87]. Bi2Te2Se has a large bulk 
resistivity because its chemical characteristic makes the for- 
mation of crystal defects more difficult. Motivated by this 
observation, there appears other tetradymite-like M2X2Y (M 
= Bi or Sb; X and F = S, Se or Te) compounds [88-90] and 
Bi2-.tSbvTe3_ySey alloy [91-96]. The alloy system with a se- 
ries of special combinations of x and y can further improve 
the bulk resistivity, and thus provide a diverse platform for 
investigating the surface transport phenomena. 
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Figure 5 Band inversion strength! of half-HeusIer compounds calculated by LDA (a) and MBJLDA (b). The materials with open symbols indicate that there 
is not experimental reports, and the lattice constant are obtained by first-principles total energy minimization. The red squares mark the topological insulator 
candidates (after applying a proper uniaxial strain), red diamonds mark topological metals, black squares mark ordinary insulators, and black diamonds mark 
ordinary metals. The topologically nontrivial phase is below the horizontal line, while the trivial phase is above horizontal line. The materials enclosed by 
circle indicate that their topological phase predicted change by using different exchange-con'elation potential. Adapted from Ref. [34]. 



It should be noted that in this review paper these 
tetradymite-Uke materials are classified into pseudo-binary 
3DTI because their crystal and electronic structures are very 
similar to that of binary Bi2Se3 family. 

3.4 3DTI in ternary compounds 

3.4.1 Half-Heusler 

The crystal structure of ternary half-Heusler compounds is 
described by space group F43m (No. 216). The chemical for- 
mula of these materials is XYZ, where X and Y are transition 
or rare earth metals and Z a heavy element. It can be regarded 
as a hybrid compound of XZ with rock-salt structure, and 
XY and YZ with the zinc -blende structure. The band struc- 
ture of half-Heusler compounds at the F point near the Fermi 
level splits into twofold-degenerated Fe, twofold-degenerated 
F7, and fourfold-degenerated Fg states. Away from F point, 
the valence bands and conduction bands are well separated 
throughout entire BZ. 

Three research groups independently predicted that un- 
der a uniaxial strain ternary half-Heusler compounds are new 
family of 3DTI [33,97,98]. Since the low-energy electronic 
structure of half-Heusler compounds is dominated at F point 
and similar to other zinc -blende semiconductors. The identi- 
fication of topologically nontrivial feature via the band inver- 
sion picture becomes possible [97,98], just as what occurs in 
HgTe. The authors in Ref. [33] also directly calculated the 
Z2 topological invariants from the bulk band structure. The 
calculated Z2 topological invariant 1;(000) definitively con- 
firmed that some of ternary half-Heusler compounds are STI. 

The band topology calculated by first-principles method 
is sensitive to the exchange-correlation potential [34,99]. 



Hence, different potentials should be used to explore the 
topological nature. The recently developed semilocal MB- 
DLDA potential is believed to be better suited for calculating 
the topological band structure [34,99]. As shown in Fig. [5] 
the topologically nontrivial and trivial phases are located be- 
low and above the horizontal line, respectively. One can see 
that the topological phases of ScAuPb and YPdBi change by 
using different exchange-correlation potentials. 

Recently, there have appeared a few experimental works 
about the electronic structures, transport properties, and topo- 
logical phenomena of ternary half-Heusler compounds [100- 
104]. Since a large number of materials with half-Heusler 
structure possess additional properties such as magnetism 
[105] and superconductivity [106], the combination of the 
predicted topological order with ferromagnetic order and/or 
superconductive order may provide an exciting platform for 
novel quantum devices. 

3.4.2 Chalcopyrite 

The crystal structure of ternary chalcopyrite compounds is 
described by the space group I42d (No. 122) with body- 
centered tetragonal structure. The chemical formula of these 
materials is ABC2, which can be regarded as a superlattice 
of two cubic zinc -blende unit cells and the A and B cations 
are ordered on the two different sites. Since the overall struc- 
tural similarity between the ternary chalcopyrites and binary 
zinc -blende analogs, the electronic structures of the former 
one are expected to closely resemble that of latter one. There- 
fore, some of the chalcopyrite compounds may have the same 
topological class with the 3DTI HgTe. 

Feng et al. [42] predicted that a large number of chalcopy- 
rite compounds can realize the topological insulating phase 
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by exploiting adiabatic continuity of band structures and di- 
rect evaluation of the Z2 topological invariants. Comparing 
with other 3DTI with cubic symmetry, a seeming advantage 
in chalcopyrite compounds is that topologically nontrivial 
state can be formed in their native states without any external 
strain. This is because the AB cation ordering and additional 
structural modifications explicitly break the cubic symmetry. 
In chalcopyrite compounds, the band inversion strength AE is 
redefined as the energy difference between the i-orbital orig- 
inated Fe states and the VBM at the F point. Figure |6] shows 
the band topology of the chalcopyrite family of the I-III-VI2 
compounds (I = Cu, Ag, Au; III = In, Tl; V = S, Se, Te), 
as well as the II-IV-V2 compounds (II = Zn, Cd, Hg; IV — 
Ge, Sn; V = As, Sb). One can see that there are a large num- 
ber of topologically nontrivial materials. More importantly, 
many chalcopyrite topological insulators have a close lattice 
matching to several mainstream semiconductors, which is es- 
sential for a smooth integration into current semiconductor 
technology. The diverse physical properties of chalcopyrite 
semiconductors [108-111] make them appealing candidates 
for novel quantum devices. 

3.4.3 Thallium-based chalcogenides 

Thallium-based chalcogenides are narrow band-gap semi- 
conductors with rhombohedral structure described by space 
group R3m (No. 166). The elemental composition of these 
materials is TIXF2 (X = Bi, Sb and 7 = 8, Se, Te), which can 



be regarded as a sequence of hexagonally close-packed lay- 
ers with the order of -T1-Y-X-Y-. In contrast to Bi2Se3 family, 
the TIXF2 are not layered compounds because each Tl or X 
atomic layer is sandwiched by two Y atomic layers and every 
two atomic layers are coupled by strong covalent bonds. As 
a consequence, the surface electronic structure is sensitive to 
surface relaxations and different terminations. 

Although the identification of topologically nontrivial 
phase in these compounds is more complex than that of 
Bi2Se3 family, several first-principles studies have predicted 
that TlSbTe2, TlSbSe2, TlBiTe2, and TlBiSe2 are 3DTI 
[112-114]. The single Dirac cone surface state of TlBiTe2, 
and TIBiSe2 have been confirmed by ARPES measurements 
[115-118]. 

3.5 3DTI in quaternary compounds 

The search for 3DTI has already been extended to quaternary 
I2-II-IV-VI4 compounds [133,134]. These materials are de- 
scribed by space group I42m (No. 121) with body-centered 
tetragonal structure. Structurally, quaternary I2-II-IV-VI4 
compounds can be viewed as the subcompounds derived from 
chalcopyrite I-III-VI2 compounds by mutating two group-Ill 
cations to one group-II and one group-IV cations. Because 
of the structural similarity, the band topology of these com- 
pounds is expected to be nontrivial in analogy to that of chal- 
copyrite 3DTI [42]. Based on first-principles calculations, 
Chen et al. [134] and Wang et al. [133] predicted that a 
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number of I2-II-IV-VI4 compounds are indeed topologically 
nontrivial materials. More importantly, the electronic proper- 
ties of these compounds can be better altered because of the 
enhanced chemical and structural degrees of freedom. For 



example, the nontrivial band-gaps can be further increased in 
Ag2HgPbSe4 (0.047 eV) and Cu2ZnGeSe4 (0.069 eV), which 
are larger than the band-gaps in strained HgTe and some of 
chalcopyrite compounds. 



Table 1: 3DTI host materials in current literatures are classified by simple elements, binary 
compounds, pseudo-binary compounds, ternary compounds, and quaternary compounds. 
Here, the tetradymite-like materials are called pseudo-binary compounds because their 
crystal and electronic structures are very similar to that of binary Bi2Se3 family. In some 
materials with zero band-gaps, proper strains should be added to realize the topological 
insulating state. Note that ternary half-Heusler and chalcopyrite 3DTI are not list one by 
one due to their large numbers, but the readers can refer to the Fig. [sjand Fig. [6] 
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4 Conclusion and outlook 

Since the potential applications in spintronic field and quan- 
tum computation, 3DTI has no doubt become topic of interest 
in the fields of condensed matter physics and material sci- 



ence. Thousands of scientific colleagues are attracted by its 
exotic quantum phenomena, pushing this field in a rapid pace. 
Similar to other stages of development in the history of con- 
densed matter physics, for example the QH eff'ect, the emer- 
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gence of high quahty samples is of crucial importance. For 
this purpose, this paper reviews the recent progress of mate- 
rials realization in 3DTI. 

Based on topological band theory [18,19], the first- 
principles calculation plays a pivotal role in the prediction 
of 3DTI host materials. We have introduced four methods 
of searching for new 3DTI, including surface state electronic 
structure, adiabatic continuity band transformation, band in- 
version picture, and Z2 topological invariants. None of these 
methods is solely used in present literatures, but two or more 
of them are usually combined together to predict topologi- 
cally nontrivial materials. To date, many materials have been 
predicted to be 3DTI (see Table [TJ. Some of them have been 
confirmed by experiments via the ARPES measurements of 
surface Dirac electronic structure. The most favorable ma- 
terials are thus the Bi2Se3 family, which have been exten- 
sively studied from bulk to thinfilm [58-71]. But remarkably 
conducting bulk state even in high quality sample blocks fur- 
ther studies of the surface transport. For decreasing the bulk 
carrier concentration, the pseudo-binary compound derived 
from Bi2Se3 family, for example Bi2Te2Se, is found to be new 
class of 3DTI with much larger bulk resistivity [86-94]. The 
theoretical prediction of 3DTI was also extended to ternary 
and quaternary compounds, but there are rare experimental 
reports about these compounds mainly because of the diffi- 
culties in their synthesis and growth. 

We propose several interesting and valuable directions in 
the further work, (i) Among the present known 3DTI host 
materials, there does not exist a unique one which can suf- 
ficiently meet various requirements in experimental obser- 
vations. Therefore, searching for new 3DTI host materi- 
als, particularly for those with simpler surface electronic 
structure and easier synthetic and growing condition, is still 
needed. This should be a central task in this rapidly devel- 
oping field, (ii) Although the present 3DTI host materials 
are often discussed in the context of the non-interacting band 
theory, the materials realization of topological Mott insula- 
tor [135,136] and topological Kondo insulator [137,138] with 
strong electron-electron interaction should also be taken into 
account. Furthermore, the topological order affected by dis- 
order, i.e., topological Anderson insulator [139-141], is in- 
teresting and needs to be thoroughly studied, (iii) The real- 
ization of the quantum anomalous Hall effect and topologi- 
cal superconductor is also a frontier direction. For quantum 
anomalous Hall effect, some theoretical models were pro- 
posed by means of magnetic doping in the film or surface 
of 3DTI [142-147], but they have not been reproduced by 
current experiments. For topological superconductor, one ex- 
pects to realize it by constructing proper interface between 
3DTI and conventional superconductors due to the proxim- 
ity effect [18,148]. In addition, a mysterious particle obeyed 
non-Abelian statistics, the Majorana fermion, could be de- 
tected in the interface and it will be hopefully utilized in the 
topological quantum computations. 



The authors were supported by National Basic Research Program of China 
(973 Program Grants No. 201 ICBAOOlOO) and National Natural Science 
Foundation of China (Grants No. 10974231 and 11174337). 

1 Wen X G. Topological orders and edge excitations in fractional quantum 
Hall states. Adv Phys, 1995, 44: 405-473 

2 Thouless D J, Kohmoto M, Nightingale M P, Nijs M d. Quantized hall 
conductance in a two-dimensional periodic potential. Phys Rev Lett, 
1982, 49: 405^08 

3 Klitzing K v, Dorda G, Pepper M. New method for high-accuracy de- 
termination of the fine-structure constant based on quantized Hall resis- 
tance. Phys Rev Lett, 1980, 45: 494^97 

4 Bernevig B A, Zhang S C. Quantum spin Hall effect. Phys Rev Lett, 
2006, 96: 106802 

5 Kane C L, Mele E J. A new spin on the insulating state. Science, 2006, 
314: 1692-1693 

6 Kane C L, Mele E J. Quantum spin Hall effect in graphene. Phys Rev 
Lett, 2005, 95: 226801 

7 Kane C L, Mele E J. Z2 topological order and the quantum spin Hall 
effect. Phys Rev Lett, 2005, 95: 146802 

8 Kohn W. Theory of the insulating state. Phys Rev, 1964, 133: A171- 
A181 

9 Yao Y G, Ye F, Qi X L, et al. Spin-orbit gap of graphene: First-principles 
calculations. Phys Rev B, 2007, 75: 041401 

10 Huertas-Hemando D, Guinea F, Brataas A. Spin-orbit coupling in curved 
graphene, fullerenes, nanotubes, and nanotube caps. Phys Rev B, 2006, 
74: 155426 

1 1 Min H, Hill J E, Sinitsyn N A, et al. Intrinsic and Rashba spin-orbit 
interactions in graphene sheets. Phys Rev B, 2006, 74: 165310 

12 Bernevig B A, Hughes T L, Zhang S C. Quantum spin Hall effect and 
topological phase transition in HgTe quantum wells. Science, 2006, 314: 
1757-1761 

13 Konig M, Wiedmann S, Briine C, et al. Quantum spin Hall insulator state 
in HgTe quantum wells. Science, 2006, 318: 766-770 

14 Fu L, Kane C L, Mele E J. Topological insulators in three dimensions. 
Phys Rev Lett, 2007, 98: 106803 

15 Moore J E, Balents L. Topological invariants of time-reversal-invariant 
band structures. Phys Rev B, 2007, 75: 121306 

16 Roy R. Topological phases and the quantum spin Hall effect in three di- 
mensions. Phys Rev B, 2009, 79: 195322 

17 Fu L, Kane C L. Topological insulators with inversion symmetry. Phys 
Rev B, 2007, 76: 045302 

18 Qi X L, Zhang S C. Topological insulators and superconductors. Rev 
Mod Phys, 2011.83: 1057-1110 

19 Qi X L, Zhang S C. The quantum spin Hall effect and topological insula- 
tors. Physics Today, 2010, 63: 33-38 

20 Hasan M Z, Kane C L. Colloquium : Topological insulators. Rev Mod 
Phys, 2010, 82: 3045-3067 

21 Hasan M Z, Moore J E. Three-dimensional topological insulators. Ann 
Rev Cond Matter Phys, 201 1, 2: 55-78 

22 Moore J E. The birth of topological insulators. Nature (London), 2010, 
464: 194-198 

23 Brumfiel G, Topological insulators: Star material. Nature (London), 
2010, 466: 310-311 

24 DiVincenzo D P, Mele E J. Self-consistent effective-mass theory for in- 
tralayer screening in graphite intercalation compounds. Phys Rev B, 
1984, 29: 1685-1694 

25 Semenoff G W. Condensed-matter simulation of a three-dimensional 
anomaly. Phys Rev Lett, 1984, 53: 2449-2452 

26 Appelbaum Joel A, Hamann D R. The electronic structure of solid sur- 
faces. Rev Mod Phys, 1976, 48: 479^96 

27 Hoffmann R. A chemical and theoretical way to look at bonding on sur- 



12 Feng W X, er at Sci Chiiia-Phys Mech Astron December (2012) Vol.55 No. 12 



faces. Rev Mod Phys, 1988, 60: 601-628 

28 Hoffmann R. Solids and surfaces: a chemist's view of bonding in ex- 
tended stractures. New York: VCH Publishers, 1988. 

29 Liu C C, Feng W X, Yao Y G. Quantum spin Hall effect in silicene and 
two-dimensional germanium. Phys Rev Lett, 2011, 107: 076802 

30 Teo J C Y, Fu L, Kane C L. Surface states and topological invariants in 
three-dimensional topological insulators: Application to Bii-^^Sb^. Phys 
Rev B, 2008, 78: 045426 

3 1 Zhang H J, Liu C X, Qi X L, et al. Electronic structures and surface states 
of the topological insulator Bii_ ,-Sb;t. Phys Rev B, 2009, 80: 085307 

32 Delin A, Kliiner T. Excitation spectra and ground-state properties from 
density-functional theory for the inverted band-structure systems y6-HgS, 
HgSe, and HgTe. Phys Rev B, 2002, 66: 035117. 

33 Xiao D, Yao Y G, Feng W X, et al. Half-Heusler compounds as a new 
class of three-dimensional topological insulators. Phys Rev Lett, 2010, 
105: 096404 

34 Feng W X, Xiao D, Zhang Y, et al. Half-Heusler topological insulators: 
A first-principles study with the Tran-Blaha modified Becke-Johnson 
density functional. Phys Rev B, 2010, 82: 235121 

35 Zhang H J, Liu C X, Qi X L, et al. Topological insulators in Bi2Se3, 
BiaTes and Sb2Te3 with a single Dirac cone on the surface. Nature Phys, 
2009, 5: 438^42 

36 Zhang X, Zhang H J, Wang J, et al. Actinide topological insulator mate- 
rials with strong interaction. Science, 2012, 335: 1464-1466 

37 Fukui T, Hatsugai Y. Quantum spin Hall effect in three dimensional mate- 
rials: Lattice computation of Z2 topological invariants and its appUcation 
to Bi and Sb. J. Phys. Soc. Japan, 2007, 76: 053702 

38 Feng W X, Wen J, Zhou J J, et al. First-principles calculation of Z2 topo- 
logical invariants within the FP-LAPW formalism. Comp Phys Comm, 
2012, 183: 1849-1859 

39 Fu L, Kane C L. Time reversal polarization and a Z2 adiabatic spin pump. 
Phys. Rev. B, 2006, 74: 195312 

40 King-Smith R D, Vanderbilt D. Theory of polarization of crystalline 
solids. Phys Rev B, 1993, 47: 1651-1654 

41 Resta R. Macroscopic polarization in crystalline dielectrics: the gemetric 
phase approach. Rev Mod Phys, 1994, 66: 899-915 

42 Feng W X, Xiao D, Ding J, et al. Three-dimensional topological insu- 
lators in I-III-VI2 and II-IV-V2 chalcopyrite semiconductors. Phys Rev 
Lett, 2011, 106: 016402 

43 Soluyanov A A, Vanderbilt D. Computing topological invariants without 
mversion symmetry. Phys Rev B, 2011,83: 235401 

44 Yu R, Qi X L, Bemevig A, et al. Equivalent expression of Z2 topologi- 
cal invariant for band insulators using the non-AbeUan Berry connection. 
PhysRevB, 2011, 84: 075119 

45 Singh D J. Planewaves, Pseudopotentials and the LAPW Method. 
Boston: Kluwer Academic, 1994. 

46 Tran F, Blaha P. Accurate band gaps of semiconductors and insulators 
with a semilocal exchange-correlation potential. Phys Rev Lett, 2009, 
102: 226401 

47 Blaha P, Schwarz K, Madsen G, et al. Wien2k, an augmented plane wave 
plus local orbitals program for calculating crystal properties. Austria: 
Vienna University of Technology, Vienna, 2001 . 

48 Groves S, Paul W. Band structure of gray tin. Phys Rev Lett, 1963, 11: 
194-196 

49 Chelikowsky J R, Cohen M L. Nonlocal pseudopotential calculations for 
the electronic structure of eleven diamond and zinc-blende semiconduc- 
tors. Phys Rev B, 1976, 14: 556-582 

50 Perdew J P, Burke K, Emzerhof M. Generalized gradient approximation 
made simple. Phys Rev Lett, 1996, 77: 3865-3868 

51 Hsieh D, Wray L, Qian D, et al. Direct observation of spin-polarized sur- 
face states in the parent compound of a topological insulator using spin- 
and angle-resolved photoemission spectroscopy in a Mott-polarimetry 
mode. New J Phys, 2010, 12: 125001 



52 Hsieh D, Xia Y, Wray L, et al. Observation of unconventional quantum 
spin textures in topological insulators. Science, 2009, 323: 919-922 

53 Lemer L S, Cuff K F, Williams L R. Energy-band parameters and relative 
band-edge motions in the Bi-Sb alloy system near the semimetallsemi- 
conductor transition. Rev Mod Phys, 1968, 40: 770-775 

54 Liu Y, Allen R E. Electronic structure of the semimetals Bi and Sb. Phys 
Rev B, 1995, 52: 1566-1577 

55 Hsieh D, Qian D, Wray L, et al. A topological Dirac insulator in a quan- 
tum spm Hall phase. Nature (London) 2008, 452: 970-974 

56 Guo H, Sugawara K, Takayama A, et al. Evolution of svuface states in 
Bii_vSbj alloys across the topological phase transition. Phys Rev B, 
2011, 83: 201104 

57 Nishide A, Taskin A A, Takeichi Y, et al. Direct mapping of the spin- 
filtered surface bands of a three-dimensional quantum spin Hall insulator. 
Phys Rev B, 2010, 81: 041309 

58 Xia Y, Qian D, Hsieh D, et al. Observation of a large-gap topological- 
insulator class with a single Dirac cone on the surface. Nature Phys. 
2009, 5: 398^02 

59 Chen Y L, Analytis J G, Chu J H, et al. Experimental realization of 
a three-dimensional topological insulator, Bi2Te3. Science 2009, 325: 

178-180 

60 Hsieh D, Xia Y, Qian D, et al. Observation of time-reversal-protected 
single-Dirac-cone topological-insulator states in Bi2Te3 and Sb2Te3. 
Phys Rev Lett, 2009, 103: 146401 

6 1 Zhang T, Cheng P, Chen X, et al. Experimental demonstration of topolog- 
ical surface states protected by time-reversal symmetry. Phys Rev Lett, 

2009, 103: 266803 

62 Chen Y L, Chu J H, Analytis J G, et al. Massive Dirac fermion on the sur- 
face of a magnetically doped topological insulator. Science 2010, 329: 
659-662 

63 Hsieh D, Xia Y, Qian D, et al. A tunable topological insulator in the spin 
hehcal Dirac transport regime. Nature (London) 2009, 460: 1101-1105 

64 Kuroda K, Arita M, Miyamoto K, et al. Hexagonally deformed Fermi 
surface of the 3D topological insulator Bi2Se3. Phys Rev Lett, 2010, 
105: 076802 

65 Zhang J L, Zhang S J, Weng H M, et al. Pressure-induced superconduc- 
tivity in topological parent compound Bi2Te3. PNAS, 2011, 108: 24-28 

66 Zhang W, Yu R, Zhang H J, et al. First-principles studies of the three- 
dimensional strong topological insulators Bi2Te3, Bi2Se3 and Sb2Te3. 
New J Phys, 2010, 12: 065013 

67 Park K, Heremans J J, Scarola V W, et al. Robustness of topologicaUy 
protected surface states in layering of Bi2Te3 thin films. Phys Rev Lett, 

2010, 105: 186801 

68 Yu R, Zhang W, Zhang H J, et al. Quantized anomalous Hall effect in 
magnetic topological insulators. Science, 2010, 329: 61-64 

69 Liu C X, Qi X L, Zhang H J, et al. Model Hamiltonian for topological 
insulators. Phys Rev B, 2010, 82: 045122 

70 Liu C X, Zhang H J, Yan B H, et al. Oscillatory crossover from two- 
dimensional to three-dimensional topological insulators. Phys Rev B, 
2010, 81: 041307 

7 1 Song J H, Jin H, Freeman A J. Interfacial Dirac cones from alternating 
topological invariant superlattice structures of Bi2Se3. Phys Rev Lett, 
2010, 105: 096403 

72 Briine C, Liu C X, Novik E G, et al. Quantum HaU effect from the topo- 
logical surface states of strained bulk HgTe. Phys Rev Lett, 2011, 106: 
126803 

73 Virot F, Hayn R, Richter M, et al. Metacinnabar (yS-HgS): A strong 3D 
topological insulator with highly anisotropic siuface states. Phys Rev 
Lett, 2011, 106: 236806 

74 Madelung O. Semiconductors: Data handbook. Berlin: Springer, 2004. 

75 Feng W X, Zhu W G, Weltering H H, et al. Strain tuning of topological 
band order in cubic semiconductors. Phys Rev B, 2012, 85: 195114 

76 Zhang W, Yu R, Feng W X, et al. Topological aspect and quantum mag- 



Feng W X, et al. 



Sci China-Phys Mech Astron December (2012) Vol.55 No. 12 



13 



netoresistance of/3-Ag2Te. Phys Rev Lett, 2011, 106: 156808 

77 Sun Y, Chen X Q, Franchini C, et al. Strain-driven onset of nontriv- 
ial topological insulating states in Zintl Sr2X compounds (X = Pb, Sn). 
Phys Rev B, 2011,84: 165127 

78 Zhu Z Y, Cheng Y C, Schwingenschlogl U. Topological phase transition 
in layered GaS and GaSe. Phys Rev Lett, 2012, 108: 266805 

79 Wang Z, Sun Y, Chen X Q, et al. Dirac semimetal and topological phase 
transitions in AjBi (A = Na, K. Rb). Phys Rev B, 2012, 85: 195320 

80 Hor Y S, Richardella A, Roushan P, et al. p-type Bi^Ses for topological 
insulator and low-temperature thermoelectric applications. Phys Rev B, 

2009, 79: 195208 

8 1 Checkelsky J G, Hor Y S, Liu M H, et al. Quantum interference in macro- 
scopic crystals of nonmetallic Bi2Se3. Phys Rev Lett, 2009, 103: 246601 

82 Analytis J G, Chu J H, Chen Y L, et al. Bulk Fermi surface coexistence 
with Dirac surface state in Bi2Se3: A comparison of photoemission and 
Shubnikov-de Haas measurements. Phys Rev B, 2010, 81: 205407 

83 Chen J, Qin H J, Yang F, et al. Gate-voltage control of chemical potential 
and weak antilocalization in Bi2Se3. Phys Rev Lett, 2010, 105: 176602 

84 Chen J, He X Y, Wu K H, et al. Tunable surface conductivity in BioSes 
revealed in diffusive electron transport. Phys Rev B, 201 1. 83: 241304 

85 Checkelsky J G, Hor Y S, Cava R J, et al. Bulk band gap and surface state 
conduction observed in voltage-tuned crystals of the topological insulator 
Bi2Se3. Phys Rev Lett, 2011, 106: 196801 

86 Ren Z, Taskin A A, Sasaki S, et al. Large bulk resistivity and surface 
quantum oscillations in the topological insulator Bi2Te2Se. Phys Rev B, 

2010, 82: 241306 

87 Xiong J, Petersen A C, Qu D X, et al. Quantum oscillations in a topolog- 
ical insulator Bi2Te2Se with large bulk resistivity (6 O. cm). Physica E, 
2012, 44: 917-920 

88 Lin H, Das T, Wray L A, et al. An isolated Dirac cone on the surface 
of ternary tetradymite-like topological insulators. New J Phys, 201 1, 13: 
095005 

89 Wang L L, Johnson D D. Ternary tetradymite compounds as topological 
insulators. Phys Rev B, 2011, 83: 241309 

90 Neupane M, Xu S Y Wray L A, et al. Topological surface states and 
Dirac point tuning in ternary topological insulators. Phys Rev B, 2012, 
85: 235406 

91 Ji H W, AUred J M, Fuccillo M K, et al. Bi2Te,.6Si.4: A topological 
insulator in the tetradymite family. Phys Rev B, 2012, 85: 201 103 

92 Zhang J S, Chang C Z, Zhang Z C, et al. Band structure engineering in 
(Bi|_vSbj)2Te3 ternary topological insulators. Nature Commun, 2011, 
2: 574 

93 Kong D S, Chen Y L, Cha J J, et al. Ambipolar field effect in the 
ternary topological insulator (Bi,Sbi_v)2Te3 by composition tuning. Na- 
ture Nanotech, 2011, 6: 705-709 

94 Arakane T, Sato T, Souma S, et al. Tunable Dirac cone in the topological 
insulator Bi2-.vSbvTe3_vSe, . Nature Commun, 2012, 3: 636 

95 Ren Z, Taskin A A, Sasaki S, et al. Optimizing Bi2-ASbvTe3-vSeY solid 
solutions to approach the intrinsic topological insulator regime. Phys Rev 
B, 2011,84: 165311 

96 Taskin A A, Ren Z, Sasaki S, et al. Observation of Dirac holes and elec- 
trons in a topological insulator. Phys Rev Lett, 201 1, 107: 016801 

97 Chadov S, Qi X L, Kiibler J, et al. Tunable multifunctional topological 
insulators in ternary Heusler compounds. Nature Mater, 2010, 9: 541- 
545 

98 Lin H, Wray L A, Xia Y, et al. Half-Heusler ternary compounds as new 
multifunctional experimental platforms for topological quantum phe- 
nomena. Nature Mater 2010, 9: 546-549 

99 Al-Sawai W, Lin H, Markiewicz R S, et al. Topological electronic 
structure in half-Heusler topological insulators. Phys Rev B, 2010, 82: 
125208 

100 Gofryk K, Kaczorowski D, Plackowski T, et al. Magnetic and transport 
properties of rare-earth-based half-Heusler phases RPdBi: Prospective 



systems for topological quantum phenomena. Phys Rev B. 2011, 84: 
035208 

101 Liu C, Lee Y, Kondo T, et al. Metallic surface electronic state in half- 
Heusler compounds RPtBi (R = Lu, Dy, Gd). Phys Rev B, 2011, 83: 
205133 

102 Butch N P, Syers P, Kirshenbaum K, et al. Superconductivity in the topo- 
logical semimetal YPtBi. Phys Rev B, 2011 84: 220504 

103 Shekhar C, Ouardi S, Fecher G H, et al. Electronic structure and lin- 
ear magnetoresistance of the gapless topological insulator PtLuSb. Appl 
Phys Lett, 2012, 100: 252109 

104 Miyawaki T, Sugimoto N, Fukatani N, et al. Structural and electrical 
propeilies of half-Heusler La-Pt-Bi thin films grown by 3-source mag- 
netron co-sputtering. J Appl Phys, 2012, 111: 07E327 

105 Canfield P C, Thompson J D, Beyermann W P, et al. Magnetism and 
heavy fermion-like behavior in the RBiPt series. J Appl Phys, 1991, 70: 
5800-5802 

106 Goll G, Marz M, Hamann A, et al. Thermodynamic and transport prop- 
erties of the non-centrosymmetric superconductor LaBiPt. Physica B, 
2008, 403: 1065-1067 

107 Lin H, Wray L A. Xia Y, et al. Single-Dirac-cone Z2 topological insulator 
phases in distorted Li2AgSb-class and related quantum critical Li-based 
spin-orbit compounds. E-print at arXiv: 1004.0999 

108 Shay J L, Wemick J H. Ternary chalcopyrite semiconductors: Growth, 
electronic properties and applications. Oxford: Pergamon Press, 1975. 

109 Medvedkin G A, Ishibashi T, Nishi T, et al. Room temperature feiTomag- 
netism in novel diluted magnetic semiconductor (Cdi_vMni)GeP2. Jpn 
J Appl Phys, 2000, 39: L949-L951 

110 Cho S, Choi S, Cha G B, et al. Room-temperature ferromagnetism in 
(Zni__vMn.,-)GeP2 semiconductors. Phys Rev Lett, 2002, 88: 257203 

111 Erwin S C, Zutic L Tailoring ferromagnetic chalcopyrites. Nature Mater, 
2004, 3: 410^14 

112 Lin H, Markiewicz R S, Wray L A, et al. Single-Dirac-cone topological 
surface states in the TlBiSe2 class of topological semiconductors. Phys 
Rev Lett, 2010, 105: 036404 

113 Yan B H, Liu C X, Zhang H J, et al. Theoretical prediction of topological 
insulators in thallium-based III-V-VI2 ternary chalcogenides. Europhys 
Lett, 2010, 90: 37002 

114 Eremeev S V, Bihlmayer G, Vergniory M, et al. Ab initio electronic 
structure of thallium-based topological insulators. Phys Rev B, 2011, 
83: 205129 

115 Kuroda K, Ye M, Kimura A, et al. Experimental realization of a three- 
dimensional topological insulator phase in ternary chalcogenide TlBiSe2. 
Phys Rev Lett, 2010, 105: 146801 

116 Sato T, Segawa K, Guo H, et al. Direct evidence for the Dirac-cone topo- 
logical surface states in the ternary chalcogenide TlBiSe2. Phys Rev Lett, 
2010, 105: 136802 

117 Xu S Y, Xia Y, Wray L A, et al. Topological phase transition and texture 
inversion in a tunable topological insulator. Science, 201 1, 332: 560-564 

118 Chen Y L, Liu Z K, Analytis J G, et al. Single Dirac cone topological 
suii'ace state and unusual thermoelectric property of compounds from a 
new topological insulator family. Phys Rev Lett, 2010, 105: 266401 

119 Yan B H. Zhang H J. Liu C X, et al. Theoretical prediction of topologi- 
cal insulator in ternary rare earth chalcogenides. Phys Rev B, 2010, 82: 
161108 

120 Yan B H, Miichler L, Qi X L, et al. Topological insulators in filled skut- 
terudites. Phys Rev B, 2012, 85: 165125 

121 Sun Y Chen X Q, Yunoki S, et al. New family of three-dimensional 
topological insulators with antiperovskite structure. Phys Rev Lett, 2010, 
105: 216406 

122 Yang K, Setyawan W, Wang S D, et al. A search model for topological 
insulators with high-throughput robustness descriptors. Nature Mater, 
2012, 11: 614-619 

123 Zhang H J, Chadov S, Miichler L, et al. Topological insulators in ternary 



14 FengWX.erai. Sci China-Phys Mech Astron December (2012) Vol.55 No. 12 



compounds with a honeycomb lattice. Phys Rev Lett, 201 1, 106: 156402 

124 Bahramy M S, Yang B J. Arita R, et al. Emergence of non- 
centrosymmetric topological insulating phase in BiTel under pressure. 
Nature Commun, 2012, 3: 679 

125 Kim J, Kim J, Jhi S H. Prediction of topological insulating behavior in 
crystalline Ge-Sb-Te. Phys Rev B, 2010, 82: 201312 

126 Sa B S, Zhou J, Song Z T, et al. Pressure-induced topological insulating 
behavior in the ternary chalcogenide Ge2Sb2Te5. Phys Rev B, 201 1, 84: 
085130 

127 Sa B, Zhou J, Sun Z, et al. Strain-induced topological insulating behavior 
in ternary chalcogenide Ge2Sb2Te5. Europhys Lett, 2011, 97: 27003 

128 Jin H, Song J H, Freeman A J, et al. Candidates for topological insula- 
tors: Pb-based chalcogenide series. Phys Rev B, 2011, 83: 041202 

129 Kuroda K, Miyahara H, Ye M, et al. Experimental verification of 
PbBi2Te4 as a 3D topological insulator. Phys Rev Lett, 2012, 108: 
206803 

130 Souma S, Eto K, Nomura M, et al. Topological surface states in lead- 
based ternary telluride Pb(Bii_.,-Sbv)2Te4. Phys Rev Lett, 2012, 108: 
116801 

131 Eremeev S V, Landolt G, Menshchikova T V, et al. Atom-specific spin 
mapping and buried topological states in a homologous series of topo- 
logical insulators. Nature Commun, 2012, 3: 635 

132 Xu S Y, Wray L A, Xia Y, et al. Discovery of several large fami- 
lies of topological insulator classes with backscattering-suppressed spin- 
polarized single-dirac-cone on the surface. E-print at arXiv: 1007.5 1 1 1 

133 Wang Y J, Lin H, Das T, et al. Topological insulators in the quater- 
nary chalcogenide compounds and ternary famatinite compounds. New 
J Phys, 2011, 13: 085017 

134 Chen S Y Gong X G, Duan C G, et al. Band structure engineering of 
multinary chalcogenide topological insulators. Phys Rev B, 2011, 83: 
245202 

135 Raghu S, Qi X L, Honerkamp C, et al. Topological Mott insulators. Phys 



Rev Lett, 2008, 100: 156401 

136 Zhang Y, Ran Y Vishwanath A. Topological insulators in three dimen- 
sions from spontaneous symmetry breaking. Phys Rev B, 2009, 79: 
245331 

137 Dzero M, Sun K, Galitski V, et al. Topological Kondo insulators. Phys 
Rev Lett, 2010, 104: 106408 

138 Dzero M, Sun K, Coleman P, et al. Theory of topological Kondo insula- 
tors. Phys Rev B, 2012, 85: 045130 

139 Li J, Chu R L, Jain J K, et al. Topological Anderson insulator. Phys Rev 
Lett, 2009, 102: 136806 

140 Groth C W, Wimmer M, Akhmerov A R, et al. Theory of the topological 
Anderson insulator. Phys Rev Lett, 2009, 103: 196805 

141 Guo H M, Rosenberg G, Refael G, et al. Topological Anderson insulator 
in three dimensions. Phys Rev Lett, 2010, 105: 216601 

142 Haldane EDM. Model for a quantum Hall effect without landau levels: 
Condensed-matter realization of the "parity anomaly". Phys Rev Lett, 
1988,61: 2015-2018 

143 Liu C X, Qi X L, Dai X, et al. Quantum anomalous Hall effect in 
Hgi_vMn,.Te quantum wells. Phys Rev Lett, 2008, 101: 146802 

144 Yu R, Zhang W, Zhang H J, et al. Quantized anomalous Hall effect in 
magnetic topological insulators. Science, 2010, 329: 61-64 

145 Qiao Z H, Yang S Y A, Feng W X, et al. Quantum anomalous Hall effect 
in graphene from Rashba and exchange effects. Phys Rev B, 2010, 82: 
161414 

146 Ding J, Qiao Z H, Feng W X, et al. Engineering quantum anoma- 
lous/valley Hall states in graphene via metal-atom adsorption: An ab- 
initio study. Phys Rev B. 201 1, 84: 195444 

147 Zhang H J, Zhang X, Zhang S C. Quantum Anomalous Hall Effect in 
Magnetic Topological Insulator GdBiTes. E-print at ai'Xiv: 1108.4857 

148 Fu L, Kane C L. Superconducting proximity effect and Majorana 
fermions at the surface of a topological insulator. Phys Rev Lett, 2008, 
100: 096407 



